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Abstract 

We present a modified TREESPH code to model galaxies in 3d. The model includes a multi-phase 
description of the interstellar medium which combines two numerical techniques. A diffuse warm/hot 
gas phase is modelled by SPH while a sticky particle scheme is used to represent a cloudy medium. 
Interaction processes, such as star formation and feedback, cooling and mixing by condensation and 
evaporation, are taken into account. Here we apply our model to the evolution of a Milky Way 
type galaxy. After an initial stage, a quasi-equilibrium state is reached. It is characterised by a star 
formation rate of ~1 Mq yr~^. Condensation and evaporation rates are in balance at 0.1-1 Mq yr~^. 

Keywords: galaxies: ISM, evolution; ISM: evolution; kinematics and dynamics; methods: 
N-body simulations 



1 Introduction 

In most 3d models of galaxies the interstellar 
medium (ISM) is described as a single phase: ei- 
ther diffuse gas e.g. using smoothed parti- 
cle hydrodynamics (SPH) (Lucy, 1977; Gingold 
&; Monaghan, 1977) or as a clumpy phase e.g. 
by sticky particles as in Theis & Hensler (1993). 
Single-phase models have been successfully ap- 
plied in the context of cosmological simulations 
(Steinmetz & Miiller, 1994, 1995; Navarro, Frenk 
&: White, 1997) as well as in simulations of iso- 
lated galaxies (Friedh & Benz, 1995; Raiteri, Vil- 
lata & Navarro, 1996; Berczik, 1999). A problem 
in these simulation can be the so-called overcool- 
ing (White & Frenk, 1991; Cole et al., 1994) lead- 
ing to an overproduction of dwarf galaxies. More 
recent models use a multi-phase ISM. This allows 
a more realistic description of star formation (SF) 
and feedback processes which could be a solution 
to this problem. 

A multi-phase ISM can be accomplished by 
modifying the SPH-algorithm (Hultman & Pha- 
rasyn, 1999; Ritchie & Thomas, 2001; Springel 



& Hernquist, 2003). We follow a different ap- 
proach and combine both treatments of a one- 
phase model in a particle based code: the hot 
diffuse gas phase is described by a SPH formal- 
ism, whereas the cold molecular clouds are rep- 
resented by sticky particles. In this model, the 
two gas phases are dynamical independent. The 
coupling between the gaseous phases is achieved 
by condensation/evaporation (C/E) and by a drag 
force due to ram pressure. Energy is dissipated by 
cloud-cloud-collisions or by radiative cooling. Fur- 
thermore, stars are formed in clouds and the stars 
return mass and energy to the gas by feedback 
processes (supernovae (SNe), planetary nebulae 
(PNe)). A similar concept is followed by Semelin 
& Combes (2002) and Berczik et al. (2003). 

In SPH codes the star formation (SF) is usu- 
ally based on the Schmidt law, i.e. the SF rate 
depends on gas density to some power and a char- 
acteristic time scale. For sticky particles the SF 
can be coupled to cloud-cloud collisions, or a sin- 
gle cloud forms stars with a constant SF efficiency. 
Here, a SF scheme using a different approach is 
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Figure 1: The left diagram shows the SF scheme: (to) the cloud is inactive (no SF); (ti) an 
embedded star cluster has been formed with a local SF efficiency; (t2) the cloud is fragmented 
by SNe energy input. In the right plot the local SF efficiency depending on the mass of the star 
forming cloud and the local pressure is shown. It is based on a model proposed by Elmegreen 
& Efremov (1997). 



presented: the process of SF is treated individu- 
ally for each cloud. The SF efficiency depends on 
local properties of the ISM and the star forming 
clouds, thereby enabling self-regulation of SF by 
feedback. 

The code is described in more detail in Sec. 2. 
In Sec. 3 a model of a Milky Way type galaxy is 
presented. A short summary and future prospects 
are given in Sec. 4. 



2 Numerical Treatment 

We use a particle code to model the evolution of 
galaxies. In our code, different components of 
a galaxy are distinguished by means of particle 
types. These particle types arc stars, clouds, (dif- 
fuse) gas and dark matter particles. The gravi- 
tational force for all particles is calculated with 
a new, very fast TREE-algorithm proposed by 
Dehnen (2002). An additional external potential 
can act as a static dark matter halo. A standard 
SPH formalism is used to compute the hydrody- 
namics of the diffuse gas. A description of the 
SPH method can be found in Monaghan (1992). 
Furthermore, the code includes the sticky particle 
scheme of Theis Sz Hensler (1993) to describe the 
collisional cloud dynamics. The integration of the 
system is done with a leap-frog scheme. Individual 
time steps are used for each particle. 



The drag force and the mass exchange rates 
for C/E are calculated according to Cowie et 
al. (1981). We apply a mass-radius-relation for 
the clouds based on observations (e.g. Rivolo Sz 
Solomon, 1988). Effects of both processes are in- 
dividually determined for each cloud using the lo- 
cal density, temperature and velocity of the hot 
gas. Energy can be dissipated by inelastic cloud- 
cloud collisions (Theis & Hensler, 1993) or by ra- 
diative cooling. The cooling function is taken from 
Bohringer Sz Hensler (1989). 

In Fig. 1 the process of SF is drawn schemat- 
ically. The general idea is that stars are formed 
in clouds and the cloud is fragmented in the pro- 
cess. Each cloud (or fragment from a previous SF 
process) is inactive for a given time Tia (to). The 
time Tia represents a global SF time scale which 
is typically of the order of a few 100 Myr. The 
time scale of the SF process itself is neglected as 
recent observations suggest that SF occurs on the 
much faster local crossing time scale (Elmegreen, 
2000). Because not all the gas in clouds is dense 
or unstable enough for SF, Tia is the time scale 
on which clouds provide the matter that is able to 
form stars. 

Once the SF process has started, an embedded 
star cluster is formed (ti) with the SF efficiency 
e being a function of cloud mass and gas pressure 
according to Elmegreen Sz Efremov (1997). The 
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Table 1: Properties of the initial disk galaxy 
model 



component 


mass"*^ 


^particles 


disk 


0.343 


20000 


- stars 


0.274 


16000 


- clouds 


0.069 


4000 


bulge 


0.167 


10000 


DM halo 


1.94 




hot gas 


0.002 


10000 


total 


0.002 


40000 



^ in units of 10^^ M, 



energy released by SNe (type II) is calculated us- 
ing a multi-power law IMF (Kroupa et al., 1993) 
with stellar lifetimes taken into account (Rateri et 
al., 1996). This energy injection disrupts the cloud 
into smaller fragments (t2) leaving a new star par- 
ticle. The time for disruption Tdis (t2— ti in Fig. 1) 
is determined from the energy input: a self-similar 
solution is used to calculate the SN shell expan- 
sion. It is assumed that the cloud disrupts when 
the shell radius equals the cloud radius. The frag- 
ments are then placed symmetrically on the shell 
with velocities corresponding to the expansion ve- 
locity at Tdig. The energy not consumed by cloud 
disruption and kinetic energy of the ejecta is re- 
turned to the hot gas phase (SPH particles). Ad- 
ditionally, the mass returned by PNe is added to 
surrounding cloud particles also taking stellar life- 
times into account. 



3 The model and results 

Initially, a galaxy similar to the Milky Way was set 

up using the galaxy building package described by 
Kuijken & Dubinski (1995). We chose their model 
A and some parameters of the initial model are 
listed in Tab. 1. A cloudy gas phase is added 
by treating every fifth particle from the stellar 
disk ctS Si cloud. The total mass in clouds is 
Mci^tot ~ 6.9 • 10^ Mq and each cloud is randomly 
assigned a time of inactivity Tja between and 
200 Myr. Finally, a slowly rotating homogenous 
gas halo with a total mass of Mhot,tot = 2 • 10® Mq 
is added. We do not follow the chemical evolution 
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Figure 2: Shown is the temporal evolution 
(solid line) and the average (dashed line) of the 
SF rate for a model of Milky Way type galaxy. 

of the galaxy, so where a metallicity is needed we 
use solar abundances. In order to treat the feed- 
back of 'old' star particles each star particle is as- 
signed an age. For the bulge stars the age is 10 Gyr 
and for disk stars 0-10 Gyr. 

Starting from this setup the galaxy model is 
evolved for 3 Gyr with all processes switched on. 
During an initial phase (0 to ~1 Gyr) the gas 
halo collapses and expands again to reach an equi- 
librium defined by the interplay of the gravita- 
tional potential, radiative cooling, heating by SNe 
and condensation/evaporation. For the following 
2 Gyr the galaxy shows a stable evolution with 
an average SF rate of ~lM0yr^^ (Fig. 2). The 
SF rate is slowly decreasing with time due to the 
consumption of the clouds while the average SF 
efficiency is roughly constant at 5 — 10%. In the 
end, the number of particles has increased by a 
factor of about 2.5. The final number of stars and 
clouds are 63 748 and 27 038, respectively. 

The mass exchange by C/E reaches a quasi- 
equilibrium after a few 100 Myr (Fig. 3, left panel). 
Before, evaporation is the dominant process due to 
the initial collapse and heating of the diffuse gas. 
In the later evolution, the condensation rate is in 
the range of 0.1-1 M0yr~^ while the evaporation 
rate occasionally drops by few orders of magni- 
tude. These drops can be correlated with a lower 
SF rate, i.e. a lower heating by SNe. 

The initially strong evaporation can also be seen 
in the comparison of the total mass in the gas 
phases (Fig. 3, right panel). The diffuse gas mass 
increases strongly in the beginning but is roughly 
constant after that. The total cloud mass is slowly 
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Figure 3: The left diagram shows the C/E rates vs. time. After an initial phase (< 300 Myr) 
both processes reach an equilibrium at roughly lM0yr^^ (300 — 1000 Myr). The lower evap- 
oration rate at ~1 Gyr and ~2 Gyr coincides with a lower SF rate. The right plot shows the 
total mass in the cloudy and the diffuse ISM. The mass of the diffuse gas is almost constant 
after an early strong increase while the cloud mass is decreasing due to SF. 



decreasing, showing the depletion by SF. A de- 
tailed look at the diffuse gas shows also that only 
a small fraction (5%) is in the halo, while most 
gas (65%) is located in the disk (i? < 20kpc and 
z < Ikpc). Some gas (30%) has even left (at least 
temporarily) the galaxy {R > 50kpc). The av- 
erage densities in halo and disk are 3 • 10~^ cm~^ 
and 1 • 10~^ cm"'^, respectively. The gas in the 
halo is hot (> 10^ K) while the gas in the disk is 
cooler (~ 10^ K) and could be described as a warm 
phase. In fact, it could form clouds from thermal 
instabilities, so that this process should also be 
taken into account in future calculations. In the 
Milky Way the total amount of warm and cold gas 
observed in the disk is ~ 6 • 10^ M0 which is about 
a factor of 2 more than what wc sec at the end of 
our simulation. In order to get a better agreement 
we could increase the initial gas mass fraction (i.e. 
select more clouds) or add an extended gas disk 
as a reservoir for a continuous gas infall. 

The disk is stable over the full 3 Gyr of evolu- 
tion. This can be seen in Fig. 4, where a com- 
parison of the stellar surface density for the full 
model (i.e. including SF and all other processes) 
with a pure N-body model is shown (left panel). 
In the full model the surface density in the cen- 
ter (R < 2kpc) has increased significantly. The 
explanation for this is that due to dissipation 
clouds have fallen towards the center before form- 
ing stars. The surface density profiles of the clouds 



confirm this idea. The velocity dispersion profiles 
(Fig. 4, right panel) show that the disk is heating 
up a little. A comparison with the pure N-body 
model shows that some of the heating is due to 
two-body relaxation but a significant amount is 
also due to the kinetic feedback of SF. Simula- 
tions with higher initial partical numbers should 
be performed to minimize effects of two-body re- 
laxation in order to determine the heating rate due 
to feedback. 



4 Summary and future Prospects 

We presented a particle code to model the evo- 
lution of galaxies with a multi-phase description 
of the ISM including a new approach to SF. A 
model of a Milky Way type galaxy shows a rea- 
sonable behavior in terms of SF rate and stability 
of the disk. 

In the near future more simulations are needed 
to determine how the results depend on initial con- 
ditions, e.g. the initial set-up for the hot gas or 
initial particle numbers and how other parameters 
in the model, e.g. the energy feedback by SNe, af- 
fect the evolution of the model galaxy. This should 
result in a reference model, that could be used 
in simulations of interacting galaxies, e.g. to shed 
more light on trigger mechanisms of star bursts. 
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Figure 4: In the left plot the surface density of the stellar disk is plotted for different times of 
3 Gyr of evolution. The upper group of curves is for the presented model including SF and all 
other processes. For a comparison the lower group of curves (scaled by a factor of 0.1 to avoid 
overplotting) is for the pure N-body evolution of the same galaxy model. The full model shows 
an increased surface density in the center as well as in the outer parts. The right diagram shows 
from top to bottom the velocity dispersions ctr, (scaled with 0.2) and (scaled with 0.1). 
A heating of the disk can be seen. 
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